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ABSTRACT 



Context. Convergent point (CP) search methods are important tools for studying the kinematic properties of open clusters and young 
associations whose members share the same spatial motion. 

Aims. We present a new CP search strategy based on proper motion data. We test the new algorithm on synthetic data and compare it 
with previous versions of the CP search method. As an illustration and validation of the new method we also present an application to 
the Hyades open cluster and a comparison with independent results. 

Methods. The new algorithm rests on the idea of representing the stellar proper motions by great circles over the celestial sphere 
and visualizing their intersections as the CP of the moving group. The new strategy combines a maximum-likelihood analysis for 
simultaneously determining the CP and selecting the most likely group members and a minimization procedure that returns a refined 
CP position and its uncertainties. The method allows one to correct for internal motions within the group and takes into account that 
the stars in the group lie at different distances. 

Results. Based on Monte Carlo simulations, we find that the new CP search method in many cases returns a more precise solution 
than its previous versions. The new method is able to find and eliminate more field stars in the sample and is not biased towards distant 
stars. The CP solution for the Hyades open cluster is in excellent agreement with previous determinations. 

Key words, astrometry - open clusters and associations: general - open clusters and associations: individual: Hyades - stars: kine- 
matics - methods: analytical - methods: statistical 



1. Introduction 

Ever since their discovery, the existence of stellar groups with 
common space motion in the solar neighborhood has been an in- 
triguing issue whose understanding is still far from complete. 
The origin and evolution of these comoving groups of stars, 
usually referred to simply as moving groups, is explained by 
different scenarios including cluster disruption, dynamic a l ef- 
fects, and accretion even ts (lEggenl Il996t iDehnenl l2QO0t iFuxl 
l200UlNavarro et al]|2004j) . Moving groups, which are observed 
as overdensities in the velocity space and exhibit a l ow inter- 
nal v elocity dispersion, typically a few km/s or less (Mathieu 
1986), allow stu dy of the large-sc ale structure and dynamics of 
the Milky Way (lAntoia et al.l 120081) . Because of perspective ef- 
fects, the proper motions of comoving star£] appear to converge 
to a single point in the celestial sphere referred to as the conver- 
gent point (hereafter CP) of the moving cluster. The CP is im- 
portant not only for determining which stars are actual members 
of a moving cluster, but also for deriving individual kinematic 
distances of moving group members, provided that their radial 
velocities are known. This is very valuable when the trigono- 
metric parallax from the ground is not accessible and Hipparcos 



1 In fact, the proper motions of a group of stars converge to a ver- 
tex either if their space motions are parallel or if they are expanding 
uniformly from a moving point. These two dynamical states are strictly 
equivalent as far as proper motions are concerned and radial ve locity 
information is needed to distinguish between them (Blaauw 1964). 



paralla xes are not available (for recent applicatio ns of this strat- 
egy see lMamaiekH2005UBertout & Genoval2006l) . 

The first algorithms implementing a method for calculating 
the CP coord inates come from Charlier and Bohlin in 1916 (see 
ISrr^fl968h . Each of them derived an equation independently 
using the position angle of stars to determine the CP position on 
the celestial sphere. However, the constant least-square coeffi- 
cients involved in solving these equations were subject to mea- 
surement errors, thus l eading t o systematic errors in the derived 
CP co ordinates. Later. ISearel d 19451) . iPetrid d 19491) and lRomanl 
(1949) proposed different ways of correcting the Charlier and 
Bohlin equations. Their strategies represent, in a first approx- 
imation, different ways of using the position angle of stars to 
determine the CP solution. 

Another approach was proposed by iBrownl d!950t) . He in- 
troduced a reference frame with the origin approximately at the 
center of the cluster and a fundamental plane defined by drawing 
a great circle through the origin in the direction of the average 
proper motion of the cluster. The proper motion of each star is 
then resolved into two components, one parallel to the reference 
plane and the other one perpendicular to the same plane, with 
the former expected to be much larger than the latter. The coor- 
dinates of the CP were then derived by applying the method of 
maximum likelihood. 

Based on this method. fjonesl (Il97ll hereafter J71) presented 
a twofold algorithm dedicated to simultaneously selecting the 
moving group members and calculating the CP pos ition. This 
method was later improved and reformulated by Ide Br uiine 
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(1999, hereafter B99) to make full use of the Hipparcos data 
and allow for internal motions within the moving group. 

There are many methods of finding the CP coordinates of 
moving groups. Several of them use, as J71 and B99 do, the 
observed proper motions of stars and differ in the details of 
the search strateg y (see for example Makarov & RobichonliOQll 
Mak arovll2007 aj|b|). Other methods use different observed quan- 
tities, such as parallaxes and radial velocities dChen et al.|[ l997 ; 
Chereul et al. 1999; Hoogerwerf & Aguilar! [19991; lAsiain et all 
1999). The choice of which method to use depends essentially 
on the observational information that is available. The new CP 
search method that we propose in this article uses the observed 
proper motions to select the stars that belong to the moving 
group and determine the CP position. It builds on the works of 
J71 and B99, so we give in Sect. 2 some details on these methods 
that will be useful in presenting our own work. 

In the rest of this paper we present and test a new version 
of the CP search method. Section 3 describes our new algorithm 
while Sect. 4 deals with the construction of synthetic data from 
moving group simulations dedicated to testing and investigat- 
ing the performance of our algorithm in comparison with pre- 
vious ones. Section 5 presents an application of our method to 
the Hyades open cluster using Hipparcos data, and finally Sect. 6 
summarizes the results of this work. 

2. The classic CP search method 

2.1. Jones' CP search method 

In J71, Jones proposed a maximum-likelihood method for si- 
multaneously determining the CP position and the members 
of the moving group. To begin with, we consider a moving 
group of stars with known positions (a, 6) and proper motions 
(p a cos 6, fig), which are the two components of the proper mo- 
tion vector fi. The tangential velocity Vt [km/s] is written in 
terms of the parallax n [mas] of the star and its proper motion 
p. [mas/yr] as 



IV, 



A \p] 



(1) 



where A = 4.74047 km yr/s is the ratio of one astronomical unit 
in km to the number of seconds in one Julian year. Furthermore, 
the radial velocity V r [km/s] of a given star is given as a function 
of the stream space velocity V [km/s] by 



|V r | = |V| cos A, 



(2) 



where A, the angular distance between the CP coordinates 
(a cp , S cp ) and a given star of the moving group, is computed from 



cos A — sin 6 sin 8 cp + cos 5 cos 5 cp cos(a cp - a) . 



(3) 



We recall that A is also the angle between V and V r for parallel 
motions, which allows us to calculate the tangential and radial 
components of the space velocity. The two components of the 
proper motion //|[, directed parallel to the great circle that joins 
the star and the CP, while p ± , directed perpendicular to the same 
great circle (see Fig. 1), can be expressed as 



n 



sin 8 cos 8\lii a cos 6 
— cos 6 sin 6 1 1 /is 



where the position angle 9 of the CP is given by 
sm(a cp - a) 



tan# = 



cos 6 tan 5 cp - sin 6 cos(a cp - a) 



(4) 



(5) 




Fig. 1. Schematic representation of the proper motion compo- 
nents parallel and perpendicular to the great circle that contains 
the star S and the CP of the moving group. 



When strict convergence occurs, i.e., when the proper motion 
vector is directed exactly towards the CP, we have p ± — and 
p.\\ = p.. In practice, strict convergence is not attained because of 
measurement errors and the internal velocity dispersion within 
the comoving group of stars. 

The basic idea of Jones' method is to determine the CP of the 
moving group by comparing p ± for each star with its expectation 
value of zero. The procedure for finding the CP of a moving 
group with N stellar candidates starts by overlaying the sky with 
a grid. Each grid point is a CP candidate whose coordinates are 
denoted by (a cp ,6 cp ). At each grid point one computes the X 2 
value 



(6) 



7=1 



where t± - p±/cr ± is the error-weighted value of p ± for each 
star in the group. Assuming that t± is normally distributed, X 2 is 
distributed as x 2 with N — 2 degrees of freedom. Minimizing X 2 
is equivalent to maximizing the likelihood function of the com- 
puted t± values. Thus, the most likely CP is the grid point with 
the lowest X 2 value. However, the lowest X 2 could still occur by 
chance, so one should evaluate the probability e that X 2 is higher 
than its observed value. It is given by 



1 



r[kN-2)]J x i 



X 2 



A[N-2]-l-x 



dx . 



(7) 



If the computed probability is too low, the star with the highest 
\t±\ is rejected. The number of stars in the sample is corrected 
and one goes back to the first step. This procedure is repeated 
until e has reached a preset value (to be discussed later). When 
this is done, the grid point in the last iteration is defined as the CP 
of the moving group, and all remaining stars in the sample are 
considered to be group members. Thus, determining the CP of a 
moving group is linked to the selection of its members. Before 
doing the analysis discussed previously, it is necessary to elimi- 
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nate stars with low t ± values, because these stars are not rejected 
by the method. The rejection criterion is defined as 



^jul cos 2 S + p] 



< t„. 



(8) 



where f„„„ is the rejection threshold (see below). 

Searching for the CP on a grid of trial points defined on the 
sky was made necessary because of the limited computer power 
available in the early 70s. However, the grid-based search repre- 
sents a robust method of determining the CP position of a mov- 
ing group, because it avoids the sometimes complicated topol- 
ogy of the X 2 function. One disadvantage of this implementa- 
tion is that it only returns the coordinates of the CP and gives no 
information about their errors. 



2.2. De Bruijne's CP search method 

Unlike the Jones method, which does not consider any possible 
internal motions within the moving group, the "refurbished" CP 
method developed by B99 also handles this and takes advantage 
of the much more powerful computing facilities available today. 
Jones's version of the CP method is modified in three ways and 
extended to include the determination of individual membership 
probabilities for each moving group star as explained below. 

First, the method is adapted to take advantage of the 
Hipparcos data regarding the propagation of errors. Second, the 
definition of t± is modified to include the velocity dispersion in 
the moving group. The selection of stars with t± = does not al- 
low one to identify all moving group members, since the proper 
motions of some are not directed exactly towards the CP because 
of their velocity dispersion. De Bruijne's method assumes that 
the velocity dispersion <x,, [km/s] and the mean distance d [pc] 
of the group are known in advance. Thus, the one-dimensional 
velocity dispersion <x,„, [mas/yr] of the group, in proper motion 
units, is given by 



°~ int — 



l OOOer,, 
Ad 



and the new definition of t ± becomes 



t, = 



(9) 



(10) 



+ erf 



The definition of f„„„ is also modified to account for cr,, !f . It is 
written 



(ID 



+ erf 



Furthermore, the grid-based approach used by J71 is replaced 
by a direct minimization routine that returns the CP position, as 
well as its uncertainty. Finally, a membership probability p is 
assigned to each moving group member and takes the velocity 
dispersion into account. It is given by 



p = exp 



A 



O 2 , +cr 2 



(12) 



As discussed by B99, the rejection threshold for a given 
star depends not only on its individual membership probabil- 
ity, but also on the membership probability of all stars consid- 
ered. Including the velocity dispersion term in the definition of 



fj_ tends to raise the individual membership probabilities, thus 
allowing the method to include stars that are not members of the 
moving group. Therefore, stars with very low membership prob- 
ability, even if not rejected by the method, should be carefully 
analyzed. On the other hand, stars with high membership proba- 
bility (i.e. low t x ) should also be verified, because the selection 
of moving group members is biased towards distant stars that 
generally have smaller proper motions, thus smaller p ± compo- 
nents. 



3. A new CP search method 

The original implementation of the CP method considers the 
stellar proper motion vector in two separate steps: a directional 
decomposition of the proper motion into the components p.\\, p ± 
followed by a minimization routine based on the amplitude of 
p±. Our algorithm differs from previous ones by considering 
both direction and amplitude of the stellar proper motions at 
once in the minimization routine. As a result, we will see that 
the new method is less biased towards more distant stars. We 
also take into account the internal motions of group members by 
introducing an individual correction for each star that depends 
on its proper motion and accounts for the fact that stars within 
the group lie at different distances. Our method makes an initial 
guess of the CP position using a grid-based approach that is sim- 
ilar to J7 1 . Once an approximate CP position is found, it is then 
refined by an analytical minimization routine to return a more ac- 
curate solution. This procedure not only guarantees the solution 
convergence but also returns both the CP coordinates and their 
uncertainties. We present below the details of our algorithm. 

The apparent motion of cluster members over the ce- 
lestial sphere takes place along the arcs of great circles. 
The idea of visualizing the proper motion of a star over a 
great circle was first i ntroduced by Herschel in 1783 (see 
Tru mpler & Weaverlll953l) and has been use d more recently by 
Aba d et al.1 d2003l) and lAbad & Vieiral d2005l) . Let r = (x, y, z) = 
(cos a cos 6, sin a cos 6, sin 6) denote the position of a star with 
coordinates (a, 6) and proper motion {p a ,ps), in the usual or- 
thogonal equatorial coordinate system, in which the unit vectors 
x, y, z point respectively to the vernal equinox, the point on the 
equator with a = 90°, and the northern equatorial pole. We com- 
pute the angular velocity of the star as the time derivative of the 
position vector 





' - sin a s 




' - cos a sin 6 s 


r = p a cos 6 


cos a 




- sin a sin 5 




{ ) 




cos 6 , 



(13) 



The great circle representing the motion of a given star is the in- 
tersection of the celestial sphere with the plane defined by r and 
r (see Fig. 2). Whenever a group of stars with parallel motions 
exists, their great circles should intersect at two opposite points 
on the celestial sphere: the CP and its mirror point. We denote 
the mirror point of the CP in this work as the divergent point 
(DP hereafter) of the moving group. Assuming that (a cp ,5 cp ) 
are the coordinates of the CP, the coordinates of the DP are 
(ad p ,6dp) = {a C p + 180°,— S cp ). From this visualization of the 
proper motions over the entire sphere one sees that a search for 
the group vertex will in fact return two geometrically equivalent 
solutions. The distinction between both solutions depends on the 
direction of the proper motion vectors. 
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' - cos Ol sin (5 




' - sin a s 


p = r x r = p a cos 6 


sin a sin 6 


-Ms 


cos a 




cos 6 




{ o J 



The great circle of each moving group member contains the 
star, the CP (DP) and the direction of motion. It can be described 
by the vector p normal to the plane defined by r and r: 



(14) 



For each group member, the vector p defines the pole of its great 
circle on the celestial sphere. The pole representation of each star 
is unique and includes its position plus its motion. Whenever a 
group of comoving stars exists, the poles of the great circles of 
its members are located in a common plane, and the intersection 
of this plane with the celestial sphere forms another great circle. 
The poles of this polar great circle are the apex and antapex 
of the moving group, i.e., its CP and DP. In practice these two 
points are the intersections of the individual stellar great circles 
(see Fig. 2). They are calculated by forming the cross products 
of all poles of members between themselves, p, x pj, when i 
and j are group members. The distribution of intersection points 
between the polar great circles of a group of stars can highlight 
preferred directions of moti on and may be u sed to search for 



moving group members (see lAbad et al 



y be use 
12003b . 




Fig. 2. Schematic representation of the motion over a great circle 
for two stars S \ and S 2 of the moving group with proper motions 
pi and p2- The poles of the corresponding great circles are given 
by Pi and Pi. 



3.1. The new CP search algorithm 

The basic idea of our procedure is to determine the polar great 
circle that best interpolates the individual great circle poles of 
the moving group stars. The poles of the interpolated great cir- 
cle return the CP and DP solutions of the moving group. To de- 
rive a first approximation of the CP coordinates we follow the 
grid-based approach developed by J71, which is adapted to our 
strategy. For simplicity we describe the algorithm in the equato- 
rial coordinate system, but it is independent of the chosen coor- 
dinate system. 



Given a sample of N stars with known positions (a, 6), 
proper motions (p a cos 5, ps), and corresponding errors, we per- 
form the following steps: 



L. 



Compute the pole (a p , 5 P ) of the great circle for each star 
in the sample and the corresponding errors (<7V,<t>). The 
derived expressions are given in Appendix A. 
Define a grid of trial CPs on the plane of the sky and number 
each grid point (i = 1,2,3, ...,N gr id). 

Start at grid point i = 1 and assume that this point is the CP. 
The coordinates of this point are referred to as (a cp , 6 cp ). 
Calculate for each star the error-weighted value s of the or- 
thogonality error k 



The orthogonality error k is defined as 



cos S cp cos 6 P cos(a, 



10. 



sin 6 cp sin 6 P - 



cp 



(15) 



(16) 



where j runs over all stars in the sample (J = 1,2,3, ...,N). 
The derived expressions for k and o~ K are given in 
Appendix B. We assume that s is normally distributed with 
zero mean and unit variance (the validity of this hypothesis 
will be discussed in Sect. 4.5). The probability distribution 
for the star j with a given s is then 



1 



Pj 



V2^ 



CX P|-^ 



5. At the given grid point compute X 2 as 



x- 



(17) 



(18) 



;'=i 



6. Determine X 2 at each grid point by repeating steps 3-5. 

7. The total probability P for the set of calculated values is 
given by 



■n 

7=1 



e 2 



(19) 



and defines the likelihood function. With s distributed nor- 
mally, X 2 is distributed as x 2 with N -2 degrees of freedom. 
Maximizing the likelihood function is equivalent to minimiz- 
ing X 2 . Define the grid point with the lowest X 2 as the most 
likely CP. 

8. The lowest X 2 could still occur by chance rather than by a 
good fit between the observations and the model. As in J71, 
evaluate the probability e that X 2 will exceed the observed 
value of X 2 by chance. It is given by 



1 



e = 



IW-2)] 



Jr»oo 
X 2 



[N-2]-l-x 



e x dx , 



(20) 



where F(x) denotes the Gamma function for x > 0. 
If e < e m ,„, that is to say, if the computed value of X 2 is 
unacceptably high, reject the star with the highest |j| value, 
correct the number of stars in the sample N —* N —1, and go 
back to step 3. Otherwise continue to the next step. 
When e has reached an acceptable value (to be discussed in 
Sect. 4.4), choose the grid point considered in the last itera- 
tion as the maximum-likelihood CP and all non rejected stars 
in the sample are identified as moving group members. 
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The algorithm presented above simultaneously selects mov- 
ing group members and calculates the CP position. It is based 
only on position and proper motion data. As in J7 1 , before start- 
ing the procedure, it is necessary to reject all stars with proper 
motion data that carry poor information because of measurement 
errors. 

3.2. Correction for internal velocity dispersion 

The new CP search algorithm described above selects those stars 
in the sample as moving group members whose orthogonality er- 
ror k approaches zero. When k approaches its expected value of 
zero, the plane defined by the individual poles of stars and the 
one that passes through the CP and DP solutions are orthogo- 
nal. The poles of group members are located on the same plane 
when the stars in the group have parallel motions. However, per- 
fect parallelism will not necessarily be achieved by some mov- 
ing group members. Therefore, selecting only those stars with 
5 = will not identify all group members. Consequently, a small 
amount of deviation of s due to the velocity dispersion in the 
group should be permitted. 

The procedure that we use to allow for internal motions 
within the moving group is similar to the one proposed by B99. 
The definition of s is changed to 



(21) 



y/crl + Ak 2 ' 



where Ak is an estimate of the one-dimensional velocity disper- 
sion in the group that is translated into the proper motion scatter 
(AjUq., Afj.fi) of group stars. It is given by 



Ak = 



3k 



dfx* a 



A/4 + 



8k 



Am 



(22) 



where ju* = jj a cos 6. The proper motion scatter in each com- 
ponent is estimated by using Eq. (0 and assuming that the ve- 
locity dispersion and the mean distance to the group are known 
in advance. We constructed synthetic samples of moving groups 
(as discussed in Sect. 4) with different velocity dispersions and 
concluded that the scatter observed in the proper motion compo- 
nents is consistent with this estimate given by Eq. (0 and similar 
in both components. One should note that the term Ak proposed 
in this work to take the internal motions within the group into 
account differs from the one introduced by B99, because the for- 
mer is an individual correction applied to each star. Although 
average group parameters (velocity dispersion and distance) are 
used, the partial derivatives in Eq. (1221) depend on the stellar po- 
sition and proper motion of each group member. Because we 
use the proper motions, we account for the various distances of 
group members. This point is particularly important when deal- 
ing with moving groups that occupy a large volume in space. 

Following the procedure developed by B99, we then attribute 
a membership probability to each group member, defined as 



Pj = exp 



+ Ak 



(23) 



3.3. Refining the CP coordinates 

Once we have found an approximated CP position (a cp , 6 cp ) of 
the moving group as explained above, we refine it by implement- 
ing a direct minimization routine in two dimensions that leads 



to the following non linear least square equations that must be 
solved. We have 



dX 2 



7=1 



da. 



CP 



dX 2 A ds, 

- — = -> > s j— — = 0. 



7=1 



(24) 



(25) 



We follow the Levenber g-Marquardt method to solve these 
equati ons as described in Madsen et al. (2004) and iPress et alJ 
(119921) . The model is approximated at each iteration by a first 
order Taylor series expansion, and we use the initial guess 
(a cp , 8 cp ) of the CP position to obtain the final coordinates of 
the CP iteratively by successive approximations. This method 
returns the CP solution and its uncertainty in the form of a 2x2 
covariance matrix. The stopping criterion for our routine is de- 
fined by the magnitude of the CP uncertainties in the covariance 
matrix since iterating to machine accuracy is generally unneces- 
sary. 

To simplify the following discussion, we denote the two CP 
search methods (CPSMs) discussed in Sect. 2 (J71 and B99) 
collectively by the name classic CP search method, or classic 
CPSM, and the new method described in Sect. 3 by the name 
new CP search method, or new CPSM. We emphasize that the 
definition of the X 2 function is not the same in the classic and 
new CPSMs. We used the appropriate definition of X 2 in the 
various tests of the classic and new CPSMs that we discuss be- 
low. Our minimization routine using the analytical derivatives 
8X 2 /da cp and dX 2 /dd cp differs from the one used by B99. 



4. Monte Carlo simulations of moving groups 

CPSMs handle two tasks: the search for the most likely moving 
group members and determination of the CP position. Detection 
of group members requires the moving group to be distinguished 
from the field population by its kinematic properties, whereas the 
position of the CP is influenced by several parameters. To inves- 
tigate both abilities of the CP method, we applied the classic and 
new CPSMs to synthetic samples of moving groups with differ- 
ent configurations. In our simulations we focused mainly on the 
differences and advantages of both methods. 



4.1. Construction of synthetic datasets 

We construct our synthetic data using the Galactic coordinates 
system, which is more convenient, particularly when dealing 
with Galactic rotation. That is, we define a rectangular coordi- 
nate system in which the unit vectors u, V, w point, respectively, 
to the Galactic center, to the direction of Galactic rotation, and 
to the Galactic north pole. The Galactic coordinates (I, b) of the 
stars are randomly drawn in the sky region occupied by the as- 
sumed moving group. We vary the Galactic longitude I of the 
cluster center from 0° to 360° in steps of 60° and assume a con- 
stant latitude of b = 0°. Our choice of b is motivated by the 
fact that young clusters and star-forming regions lie very close 
to the Galactic plane. The individual distances are drawn from 
a Gaussian distribution where the mean distance depends on the 
cluster simulation. The stars are distributed in a distance range 
of 30 pc projected along the line of sight. We assume three dif- 
ferent configurations of distance and vary the projected angular 
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Table 1. Mean distance, field size, and number of moving group 
and field stars in the synthetic samples for each configuration 
(A,B,C) considered in this work. 



# 


Distance (pc) 


Field Size Ixb (deg 2 ) 


Nmg 


N field 


A 


100 


20x20 


200 


1000 


B 


200 


15 x 15 


100 


500 


C 


400 


10 x 10 


50 


250 



size of the group. Table 1 summarizes the characteristics of our 
synthetic samples. 

The stellar velocity components V = (u, v, w) are randomly 
drawn from a sphere in velocity space with radius of 20km/s. 
The three velocity components and position are translated into 
proper motion (/// cos b,fj.t>) and radial velocity V r for each star, 
after which observational errors are added. Position and proper 
motion errors are also drawn from Gaussian distributions with 
means of 1 mas and 1 mas/yr, respectively. We then add to the 
streaming motion of the group three components representing 
the velocity dispersion, the reflex of the Solar motion with re- 
spect to the local standard of rest, and the Galactic rotation, in the 
following manner. Each of the dispersion velocity components 
(6u, 6v, 5w) is independently drawn from a Gaussian distribution 
with an isotropic velocity dispersion of cr v = 1 km/s. Our choice 
of cr,, is motivated by the fact that young stellar groups often ex- 
hibit a low-velocity dispersion. For th e solar motion we adopt 
(U,V,W) Q = (10.00,5.23,7.17) km/s dDehnen & Binnevll 19981) 
B and A = 14.82 km s^kpc" 1 and B = -12 37 kms^kpcT 1 
for the Oort constants (Feast & Whitelocklll997h . We introduce 
Galactic rotation in the simulate d stellar prop er motions by using 
the first-order formulae given in iGreenl (119851) . 

Our synthetic samples of stars include both the moving 
group stars and a field population. The density of stars in each 
field corresponds approximately to that of the Hipparcos cata- 
log (~3 stars/deg 2 ). Field stars are located in the same field of 
view and have the same observational errors as moving groups 
stars. They are unrelated to the moving group and have random 
motions in amplitude and direction. 

Each run consists of 100 Monte Carlo cluster simulations 
for each Galactic position and distance of the cluster. We follow 
the procedure adopted by B99 to correct for internal motions by 
taking the value for cr v that was actually used to construct our 
synthetic data samples. For both CPSMs, we assume the stop 
value of the probability e defined in Eq. [7] to be e mi -„ = 0.954, 
as suggested by B99, (we come back to this in Sect. 4.4). In the 
following, we investigate and compare both CPSMs using the 
synthetic data sets. We first investigate the ability of each CPSM 
to find group members and reject field stars (see Sect. 4.2), then 
we come to a more detailed analysis of the CP itself as a function 
of several parameters (see Sect. 4.3). 

4.2. Selection of cluster members and rejection of field stars 

The CP method selects moving group members based on their 
proper motions and reject field stars that are not in the same 
proper motion range. A distinction must be made between field 
stars at the same distance as the moving group and background 
stars. To begin with, we construct samples as described in 

2 Although a mo r e rece nt value of the solar motion given 
by Schonrich et al. (2010) exists, we adopted the value of 
Dehnen & Binnev ( 1998), which has been widely used in the lit- 
erature. The results and conclusions of this paper do not depend on the 
specific value used for the solar motion. 



Sect. 4.1 and consider the case where the field population is 
at the same average distance as the moving group. The results 
of these simulations are shown in Fig. 3. On the whole, both 
CPSMs exhibit a similar performance even if the fraction of 
rejected field stars by the new CPSM is higher in some spe- 
cific cases . Our simulations show that at the shorter distances 
of Table 1 the fraction of cluster members that can be retrieved 
for both CPSMs is higher than ~ 80% and the contamination of 
field stars amounts to ~ 20% of the total number of field stars. At 
greater distances, the ability to find and eliminate field stars de- 
creases for both CPSMs, as expected, since the quality of proper 
motions also diminishes. 

We now consider the case where the field population consists 
of background stars located two times farther than the moving 
group. We perform the same simulations as before and present 
the results in Fig. 4. We find that the performance of the two 
methods differ; the new CPSM is likely to eliminate more back- 
ground stars than the classic CPSM. The difference between both 
CPSMs is more evident at greater distances (e.g. at d = 400 pc) 
where the distributions of rejected field stars for each CPSM 
are clearly separated. As already discussed by B99, the classic 
CPSM is biased toward distant stars. The definition of the X 2 
function to be minimized in the classic CPSM considers only 
one directional component of the proper motion vector, which 
leads to a biased selection of stars with small fi ± components. 
However, not all stars with low yu ± are necessarily members, be- 
cause stars at farther distances, such as background stars, gen- 
erally have smaller proper motions (i.e. small fi ± ). On the other 
hand, the definition of the X 2 function using great circle poles 
in the new CPSM considers both direction and amplitude of 
the proper motion vector at the same time. By considering the 
full proper motion vector instead of only one directional compo- 
nent (as in the classic CPSM) in the definition of the minimizing 
function X 2 leads to a membership selection that is not biased 
towards distant stars. 

4.3. Convergent point analysis 

As discussed in Sect. 4.2 the CP method is able to identify the 
moving group with a contamination of field stars that amounts 
to ~ 20% of the total field population in the sample. Although 
the method retrieves a high fraction of cluster members, the CP 
is shifted from its true position due to the contribution of the re- 
maining field stars in the sample. To better investigate the CP 
solution itself it is necessary to work on a sample where the 
moving group is the dominating population. In the following, 
we set Nmg - 80 and Nfutd - 20 in all configurations (A,B,C) 
of Table 1, construct our samples as described in Sect. 4.1, and 
present the results of CP analysis. 

4.3.1 . The effect of velocity dispersion on CP positions 

As explained earlier, the velocity dispersion of the cluster pre- 
vents the proper motion vectors from pointing exactly towards 
the CP. Here we investigate how internal velocity dispersion 
affects both CPSMs by comparing the CP positional accuracy 
found for moving groups with different values of velocity dis- 
persion. To do so, we construct three sets of synthetic data as 
described above, with cr v =1.0, 2.0, and 3.0 km/s. We display in 
Fig. 5 the CP positions obtained with both CPSMs, overlaid on 
contour lines of constant X 2 values for a cluster with no internal 
velocity dispersion, centered on the same position and distance 
as the synthetic groups with nonzero internal dispersion veloc- 
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Fig. 3. Fraction of rejected field stars for both CPSMs in the 1 800 
MC simulations. The field population is assumed to be at the 
same distance as the moving group. 



Fig. 4. Fraction of rejected field stars for both CPSMs. The re- 
sults summarize 1800 MC simulations. The field population con- 
sists of background stars at twice the distance of the moving 
group. 
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ities. The results shown in Fig. 5 are for a moving group with 
(/,*) = (180°, 0°) and d =100 pc. We find that the CP posi- 
tions derived by both CPSMs are distributed along an elongated 
structure. The great circle on which the CP lies is defined by the 
proper motion vectors of group members, however the precise 
location of the CP over this great circle appears uncertain. As 
seen in Sect. 3, when strict convergence to the CP occurs, the 
individual poles of group members lie on a common plane that 
defines the polar great circle in the celestial sphere. Because of 
internal motions within the group, the poles are scattered around 
the polar great circle, and the precise position of the CP is not 
well defined, w hich explains th e elongated structure observed in 
Fig. 5 (see also lBertiaulll958h . When comparing both CPSMs, 
we find that they perform similarly and that the dispersion of CP 
positions always grows with the increasing velocity dispersion 
within the group. 

4.3.2. Precision of the CP position derived by each CPSM 

We now discuss the CP position uncertainties as derived by the 
covariance matrix. We compare the CP precision of both CPSMs 



1 .0 km/s 



in terms of the quantity Acr cp = cr\ 



.classic 
cp 



- CT 



cp 



where cr c 



the combined uncertainty of the CP position (l cp ,b cp ) for each 
method, is given by 



cr, 



cp 



4 



LT, + CT 



b cp 



(26) 



The results of 1800 MC simulations for six values of the Galactic 
longitude (see above) and three different distances are shown in 
Fig. 6. When A<x cp (normalized in order to represent only rel- 
ative values between both CPSMs) is positive, the CP position 
uncertainties derived by the classic CPSM are higher than the 
ones obtained with the new CPSM. We find that the new CPSM 
returns a more precise CP solution for 95% of the 1800 simu- 
lations. On the other hand, we also observe that the number of 
cases where Acr cp < 0, (i.e., where uncertainties in the classic 
CPSM are lower than in the new CPSM) grows with increasing 
distance of the moving group. In other words, the new CPSM 
is less precise than the classic CPSM at large distances, where 
the higher degree of concentration of the stars affects the deter- 
mination of the CP position, defined as the intersection of the 
individual great circles of group members. 



4.3.3. Role of moving group parameters on the solution 
accuracy of the new CPSM 

The precision of the CP position is also influenced by several 
other parameters such as observational errors on proper motions, 
number of moving group members, and angular distance from 
the moving group to the CP. The most efficient way to investi- 
gate the effect of each parameter is to consider an ideal cluster 
model in which the stellar velocity components only result from 
the streaming motion of the cluster (i.e., the internal velocity dis- 
persion is zero and the effects of Galactic rotation are not consid- 
ered). We find, as expected, that for an ideal moving group both 
CPSMs return the same solution, with X 2 — and e — 1 . In the 
following, we focus on the role of each parameter for the results 
obtained with the new CPSM. We emphasize that the qualitative 
results of this investigation, obtained by using an ideal synthetic 
moving group, depend neither on the specific cluster position 
(/, 0°) nor on its distance. However, the amplitude of the CP un- 
certainties as given in Fig. 7 are specific to configuration #A of 
Table 1 with the cluster centered at (180°, 0°), which we present 
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Fig. 5. Relative position (Al cp , Ab cp ) of the CP for 100 simu- 
lations and three different velocity dispersions. Blue dots and 
red crosses denote, respectively, the results of the classic and 
new CPSMs. The contours indicate X 2 levels linearly spaced by 
AX 2 = 2000 for a cluster with the properties given in the text. 
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Fig. 6. Comparison of the CP errors between both methods in 
the 1800 MC simulations of clusters located at three different 
distances. 



here as an example to show the relative importance of each mov- 
ing group parameter in the CP solution. 

We first investigated the influence of proper motion errors on 
the CP position precision. We constructed synthetic samples of 
an ideal moving group (as defined above) by varying the relative 
error on proper motion data (in each component) from 5% to 
100%. The results are presented in Fig. 7a and show that the 
errors on the CP position grow linearly with the increasing errors 
on proper motions. 

We now come back to the distance of the moving group and 
its effect on the CP position precision. As mentioned above, we 
expect the distance to influence the CP position errors since stars 
at greater distances have smaller proper motions, and the angu- 
lar size of the group changes with respect to the distance. The 
question arises whether the precision of the CP is more affected 
by the deterioration of proper motion data at larger distances 
or by cluster concentration. Our numerical simulations make it 
possible to separate both effects. To do so, we simulate different 
moving groups varying one of these quantities while the other 
remains fixed. The results are shown in Figs. 7b and c. We find 
that cluster concentration is the dominant source of error on the 
CP position at large distances. 

We investigated the effect of the number of group members 
used for the CP position determination. Figure 7d shows the re- 
sults of this computation, where we varied the number of moving 
group members from 10 to 300. We find that above N = 200 no 
significant improvement in the CP position precision is achieved. 

Finally, we seek to determine the effect of the angular dis- 
tance from the moving group to the CP. To do so, we construct 
synthetic data sets by varying the Galactic longitude of the clus- 
ter center from 0° to 360° in steps of 0.15°. We consider the stel- 
lar three-dimensional velocity components to be isotropic, which 
fixes the CP position at (45.00°, 35.26°). Our results are shown 
in Figs. 7e and f. We find that moving groups displaced by 180° 
have the same CP errors. This is because their relative position to 
the nearest solution (CP or DP) is the same. Those groups whose 
angular distance to the CP exceed 90° are closer to the DP than 
to the CP. We recall that both solutions are always obtained with 
the same probability, so that the X 2 values are the same. We con- 
clude that the angular distance between the moving group and 
the CP influences the precision of the solution. The least precise 
CP solution is obtained for moving groups displaced by 90° with 
respect to the CP (or DP). 

4.4. Choice of ' e min 

The final selection of moving group members depends mainly 
on the choice of e m ,„, the stop parameter that should allow find- 
ing the largest number of moving group members with the least 
contamination by field stars. So far, we have adopted in our new 
CP strategy the value of e,„,„ = 0.954, which was used by B99 
in his implementation of the classic CPSM. This choice allows 
one to compare the performance of both CPSMs under the same 
conditions. To do so, we investigated whether the adopted value 
of e m i n by B99 is truly the optimum one for our new strategy. 
To do so, we ran the new CPSM on our synthetic data sets and 
varied the value of e m ,„ from 0.00 to 0.99. We then computed 
the fraction of true cluster members and field stars in the sample 
of stars selected by the new CPSM. The results are presented in 
Fig. 8 for a moving group with (l,b) = (180°, 0°) at d =100pc. 
Similar conclusions apply to all cluster configurations consid- 
ered in this work. As already expected, the relative fraction of 
group members in the sample of nonrejected stars is lower for 
the classic CPSM. This is because more field stars are accepted 
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by this method. We verified that the number of rejected group 
members is equivalent for both CP methods. We find that for 
£mi„ > 0.954 the number of moving group stars decreases while 
the fraction of field stars in the sample remains practically con- 
stant. Thus, we conclude that the value of e, m „ = 0.954 can also 
be used in the new CPSM, since it offers the best compromise 
between the number of accepted group members and rejected 
field stars. 



4.5. Normality tests applied to the s distribution 

In Sect. 3.1 we assumed that s was normally distributed with 
zero mean and unit variance. Coming back to this point, we dis- 
cuss the validity of this hypothesis. We carried out another set of 
Monte Carlo simulations with different configurations of cluster 
position, distance, angular size, and velocity dispersion in order 
to investigate whether the observed s distribution is normal with 
zero mean and unit variance. The Kolmogorov-Smirnov normal- 
ity test applied to the synthetic datasets reveals no significant 
differences between the simulations and the model. We present 
in Fig. 9 the results of this analysis for 1000 Monte Carlo simu- 
lations with different cluster configurations. We observe that the 
Kolmogorov-Smirnov test gives a p-value that is higher than the 
adopted significance level of 5% for 99.6% for our simulated 
moving groups. Thus, we conclude that the assumption of a nor- 
mal distribution of s is valid. 



5. Application to the Hyades 

Ever since the discovery that it is a nearby moving group, 
the Hyades open cluster has played a central role in astron- 
omy in studies of Galactic structure, chemical evolution, dis- 
tance calibration and stellar evolution models (see for example 
Ide Bruiine et al.ll200lL and references therein). Here we use the 
Hyades open cluster as a first application of our new CPSM 
to real data. We present the CP search using data from the 
Hipparcos catalog and compare our results with previous ones. 



5.1. Covariance matrix 

In general astrometric catalogs quote only standard errors on the 
observables. However, the Hipparcos catalog provides for each 
star a 5x5 covariance matrix for the measured astrometric pa- 
rameters (position, proper motion, and parallax) and proper use 
of the full covariance matrix is mandatory. The new CPSM is 
adapted to consider the full covariance matrix of the Hipparcos 
data, and the specific formalism used for the propagation of er- 
rors is described in Appendix C. 



5.2. CP analysis of the Hyades open cluster 

The sample of Hyades members t hat we use here consi sts of 
218 stars from Hipparcos as given in lPerrvman et al.l (fl998. here- 
after P98), whose membership analysis was based on the three- 
dimensional velocity of the stars and the structure of the cluster. 

As already discussed by B99, a nonzero velocity disper- 
sion is required in order to retrieve all cluster members. We 
adopt the value cr v =2.0 km/s used by B99 and the mean dis- 
tance of d =46 pc as given in P98. After applying the new 
CPSM to the 218 members as the starting sample, the posi- 
tion of the CP is (a cp ,6 cp ) = (97.69°, 6.30°) + (0.89°, 0.36°), 



Moving Group Members 



"D 

CD 
N 



o 



"D 

£ 

CD 
O 
O 
CO 



cz 
o 



00 

o 



CD 
d 



C\j 

d 



o 
d 



Field Stars 



Emin=0.954 , 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



Fig. 8. Fraction of accepted group members and field stars in 
the sample as a function of the adopted value for e„„„. Circles 
denote field stars and triangles cluster members. Open and 
filled symbols represent the results obtained with the classic 
and new CPSMs, respectively. The vertical dashed line marks 
em/« = 0.954. Each point represents the average value of 100 
Monte Carlo simulations for a cluster with (l,b) = (180°, 0°) 
andc/ =100pc. 
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plied to 1000 Monte Carlo simulations of moving groups with 
different configurations. The red dashed line marks the 5% level 
of significance. 
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and we retrieve 217 cluster members with a correlation coef- 
ficient of p — -0.75. Our solution is in good agreement with 
B99, (a cp ,6 cp ) = (97.81°, 6.74°) + (0.52°, 0.21°) with 213 clus- 
ter members and p = -0.84 . Another CP so lution for the Hyades 
open cluster was found bv lSchwanldl99U hereafter S91), who 
obtained (a cp ,6 cp ) = (97.68°, 5.98°) ± (0.42°, 0.18°) using only 
145 stars and proper motions from the FK5 and PPM catalogs. 
Figure 10 illustrates the proper motion vectors of the Hyades 
members and the derived CP coordinates. We find that only 
one star, HIP 28774, was rejected in our analysis as a non- 
member of the Hyades open cluster. Its proper motion is given 
as Qi a cos 6, /J S ) = (-2.46,- 7.00) + (18 .43,9.99) mas/yr in the 
original Hipparcos catalog (lESAl 119971 hereafter HIP97) and 
as (fi a cos 6, fi s ) = (42.95 -2.16)+ (14 68,8 .04) mas/yr in the 
new Hipparcos reduction (Ivan Leeuwenll2007l hereafter HIP07). 
These two values are different with low quality, as shown by 
the large uncertainties. The mean uncertainty on the proper mo- 
tions of other cluster members is < 2 mas/yr in each compo- 
nent. Furthermore, the parallax changes from7r#/p97 = (12.81 + 
12.80) mas to 7i mPQ1 = (39.18 + 11.54) mas, which reflects the 
uncertainties of the Hipparcos data for this star. That the new 
CPSM considered that HIP 28774 is not a Hyades member thus 
appears justified. 



and given in Table 2 of P98. The individual parallaxes of group 
members are given by 



Vrtan/T 



(27) 



where the corresponding quantities and units are defined in 
Sect. 2. This equation is als o valid for lin ear expansion of the 
group mentioned before (see iBlaauwil 19641) . The parallax uncer- 
tainty is derived by error propagation of this equation and takes 
the error budget of proper motions and radial velocities into ac- 
count. 

We present in Fig. 11a comparison of the parallaxes derived 
in this work with the ones given in HIP97 and HIP07. We find a 
slight bias when comparing our results with HIP97, which dis- 
appears when we compare our parallaxes with the presumably 
more accurate ones given in HIP07. Another dataset that we used 
to compare the individual parallaxes of this wo rk are the secular 
paralla xes of Hyades cluster members given bv lde Bruijne et all 
(2001). These secular parallaxes were calculated using proper 
motion data from the Hipparcos and Tycho2 catalogs. The re- 
sults of this comparison for each data set are shown in Fig. 12. 
We conclude that our results agree well with the more precise 
secular parallaxes of the Hyades open cluster. 
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Fig. 10. Proper motion vectors of the Hyades members and the 
CP coordinates derived in different works. The CP error bars are 
too small to be seen. 



5.3. Determination of individual parallaxes 

Once the CP coordinates of the Hyades have been determined, 
a final check and further extension of our CPSM lies in com- 
puting individual parallaxes and comparing it with HIP97 and 
HIP07. The determination of individual parallaxes in this work 
is restricted to group members whose radial velocity is known 



6. Summary 

We presented a new CPSM based on the idea of representing the 
stellar proper motions by great circles on the celestial sphere. 
The new CP method combines (i) maximum-likelihood analy- 
sis to simultaneously determine the CP and select moving group 
members with (ii) a direct minimization routine used to return a 
more refined CP position and its uncertainties. Our new method 
allows us to correct for internal motions within the group by ap- 
plying an individual correction for each star that depends on its 
proper motion. This procedure takes into account that the stars 
in the group are not located at the same distance. 

We performed extensive Monte Carlo simulations to test and 
compare our new CPSM with the classic CPSM regarding the 
selection of moving members and the accuracy of the CP solu- 
tion. We investigated the effects of (a) the velocity dispersion of 
the cluster, (b) observational errors on proper motion, (c) clus- 
ter distance, (d) number of group members, and (e) angular dis- 
tance to the CP on the CP solution. Our new CPSM returned a 
more precise CP solution than the classic CPSM for 95% of the 
simulations. We verified that in the absence of velocity disper- 
sion (ideal moving group) both methods exhibit the same per- 
formance. We also found that the new CPSM finds and elimi- 
nates more field stars than the classic CPSM. This situation is 
even more evident in the presence of background stars that gen- 
erally have smaller proper motions. The new CPSM is able to 
retrieve more than 80% of all cluster members with a contamina- 
tion around 20% of the total number of field stars at distances of 
100 pc and 200 pc. At larger distances the efficiency of rejecting 
field stars decreases, but it is still higher than the efficiency of the 
classic CPSM. The definition of the X 2 function in the new CP 
method considers both the amplitude and direction of the stellar 
proper motion vector instead of only one directional component 
that removes the bias towards distant stars that is observed in the 
classic CPSM. Additional information (e.g. parallaxes and radial 
velocities) is required to eliminate the remaining field stars in the 
sample that were not rejected by the CPSM. 

The new CPSM is shown to work well when applied to the 
Hyades open cluster, and its results are agree well with previous 
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Fig. 11. Individual parallaxes of Hyades members derived in this work compared with their counterparts in HIP97 (left panel) and 
HIP07 (right panel). The red line represents the expected distribution for equal results. The rms with respect to HIP97 is 3.62 mas 
and 3.34 mas for HIP07. The mean difference between the parallaxes derived in this work, and the ones in HIP97 and HIP07 are 
1.18 mas and 0.02 mas, respectively. 
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Fig. 12. Individual parallaxes of Hyades members derived in this work compared to the secular parallaxes of de Bru iine et al.l d2001l) 
derived using proper motion data from the Hipparcos (left panel) and Tycho2 (right panel) catalogs. The red line represents the 
expected distribution for equal results. The rms with respect to the Hipparcos secular parallaxes is 1.55 mas and 1.68 mas for 
Tycho2 secular parallaxes, and the mean difference between the parallaxes is 0.23 mas and -0.27 mas. 
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determinations of the CP position. The individual parallaxes de- 
rived using our CP solution are fully consistent with the trigono- 
metric parallaxes given in the Hipparcos catalog. Although our 
results agree well with data from both reductions, we find that 
our external precision is slightly better for HIP07. Our results 
are also in good agreement with the secular parallaxes found in 
the literature, which are more precise than the trigonometric par- 
allaxes given in the Hipparcos catalog. 

The original implementation of the new CPSM was adapted 
to take the full covariance matrix of the Hipparcos catalog into 
account. The same procedure can also be applied to the future 
data of the Gaia mission, which will also be published with the 
full covariance matrix. The new CPSM will be used in forth- 
coming papers to investigate the kinematic properties of several 
young associations. 
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Appendix A: Pole coordinates {a p ,5 p ) 

The rectangular coordinates of the normalized polar vector p are given as 



+ M 2 6 



- cos a sin 6 \ 

- sin a sin 6 

cos 5 



Jit 



+ M 2 6 



i - sin a '\ 
cos a 




where /v - /j a cos6. The determination of the celestial coordinates (a p ,6 p ) of the pole requires transformation to 
polar coordinate system. This yields 



a' = arctan 
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The corresponding errors (ov, ct>) are obtained by error propagation of Eqs. (lA.2b and ( IA.3b . We have to first order 
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The partial derivatives are given as 
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Appendix B: Orthogonality error 

Let us assume that e cp is a unit vector that defines the direction of the CP (a cp , 6 cp ) in the celestial sphere, and p, the normalized 
polar vector as given by Eq. dA.U . defines the pole (a p ,6 p ) of the great circle for a given star in the group. We define k = K% cp 
where k is the orthogonality error and denotes the amount of p that is projected in the e cp direction. It is written as 

k = (P • e cp ) e cp = (cos 9) e cp , (B.l) 

K 

where 9 is the angle between p and e cp . However, 9 is also the angle between the planes that contain p and e cp , which defines it as a 
spherical angle. It is given by the cosine formula in the spherical triangle that contains the CP (a cp ,d cp ) and the pole (a p , 6 P ) as 

k = cos9 = sin S cp sin 6 P + cos 6 cp cos 6 P cos(a cp - a p ). (B.2) 

Whenever perfect parallelism of the stellar motions is achieved, p and e cp are orthogonal and k = 0. In practice, one needs to 
minimize k to search for the CP and moving group members. The corresponding error of k is given by error propagation and takes 
position and proper motion errors into account. It is given as 

cr 2 . = \- cos S cp cos 6 P sin(a p - a cp )] cT 2 aP + [sin 5 cp cos 5 P - cos 6 cp sin 8 P cos(a p - a cp )] cr\„ . (B.3) 



Appendix C: Error propagation 

The Hipparcos catalog provides the five astrometric parameters (a,6,fi* a ,fj.g,n) together with the full covariance matrix. Thus, the 
propagation of errors must consider the covariances between all observables. The complete expression for the error propagation is 
given as 
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where p„ denotes the correlation coefficient of the observables x and y. The partial derivatives are given in Appendix A. 
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